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In parallel with Kibble's description of the onset of phase 
transitions in the early universe, Zurek has provided a simple 
picture for the onset of phase transitions in condensed matter 
systems, supported by agreement with experiments in 3 He 
and superconductors. In this letter we show how experiments 
with annular Josephson tunnel Junctions can, and do, provide 
further support for this scenario. 



As the early universe cooled it underwent a series of 
spontaneous phase transitions, whose potential inhomo- 
geneities (monopoles, cosmic strings, domain walls) have 
observable consequences, for structure formation in par- 
ticular. These defects appear because the correlation 
length £ of the field (or fields) whose expectation value is 
the order parameter is necessarily finite for a transition 
that is implemented in a finite time 

Using nothing more than simple causal arguments Kib- 
ble HH made estimates of this early field ordering, and 
the density of topological defects produced at GUT tran- 
sitions at 10 -35 s. Unfortunately, because the effects of 
their evolution are not visible until the decoupling of the 
radiation and matter 10 6 yrs later, it is impossible to pro- 
vide unambiguous checks of these predictions. However, 
causality is such a fundamental notion that Zurek sug- 
gested HQ that identical causal arguments, with similar 
predictions, were applicable to condensed matter systems 
for which direct experiments on defects could be per- 
formed. The hope is that successful tests of these pre- 
dictions could lead to a better understanding of phase 
transitions in quantum fields. 

Whether for the early universe or condensed matter, 
consider a quench of the system in which its temperature 
T{t) is reduced as time passes. In the vicinity of the 
critical temperature T c we assume that the temperature 
T decreases linearly with the time t at a rate dT/dt = 
—T c /tq, tq being the quenching time. 



Suppose that the 'equilibrium' correlation length 
£eg(i) = £,eq{T(t)) of the order-parameter field, and its 
relaxation time r(t), diverge at t = (when T — T c ) as 



t 



T{t) = T 



t 



(1) 



The fundamental length and time scales £0 an( i r o of a 
system are determined from its microscopic dynamics. 
One definition of r(t) is that c(t) = £ e q(t)/T(t) denotes 
the maximum speed, at time t, at which the order pa- 
rameter can change. In quantum field theory c(t) = Co, 
the speed of light in vacuo. 

Although £ e q(t) diverges at t — this is not the case for 
the true non-equilibrium correlation length £(t). Kibble 
and Zurek made two assumptions. Firstly, the correlation 
length £ of the fields that characterizes the onset of order 
is the equilibrium correlation length £ = £ e q(t) a t some 
time t constrained by causality. Secondly, at this time, 
defects appear with separation £def = 

To determine i we rephrase the original Kibble-Zurek 
argument in a way appropriate to Josephson tunnelling 
Junctions (JTJs), so as only to discuss times t > 0. 
We begin by noting that, in the adiabatic regime away 
from the transition, static defects can be thought of 
as kinks, balls, lines^] or sheets of 'false' vacuum or 
disordered ground state, of thickness 0(£ eq (t)). Thus 
£eq{t) — d£ eq (t)/dt < measures the rate at which 
these defects contract, i.e., the speed of interfaces be- 
tween ordered and disordered ground states. Since £ e q(t) 
decreases with time t > 0, the earliest possible time t 
at which defects could possibly appear is determined by 
l£eq(£)| = c(i), given our definition of c(t). Although this 
gives a lower bound for i, as an order of magnitude esti- 
mate we identify this time t with t, whence 



r « f=(r5r ) 1 /(7+i)« TQ . 
The corresponding smallest Q] correlation length is 

1-/(7+1) 
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Because of the qualitative nature of the arguments, fac- 
tors close to unity are omitted^. Further, at the same 
level of approximation, we shall use mean field critical 
indices throughout. Measurements ||bI of total vortex 
density in transitions of 3 He — B support the result (^), 
when taken together with £de/ = 0(£). 

As an independent test of the assumptions Zurek sug- 
gested using (||) to measure topological defect density (in 
which defects and anti-defects carry opposite weight). 
This is most easily done in 'one-dimensional' annular 
geometries, for which experiments were originally pro- 
posed Q with i He which, however, has £ so small that 
the creation of an effectively one-dimensional system is 
extremely difficult. A recent experiment j?Jwith annu- 
lar arrays of high-T c superconducting islands coupled by 
grain boundary Josephson junctions confirms part of the 
picture, but suffers in enforcing a predetermined domain 
structure. We also wish to check prediction (JsJ) by us- 
ing annular JTJs. As we shall see, for such JTJs £ is 
macro scopically large, permitting them to be effective 
one-dimensional systems. 

An annular JTJ consists of two superimposed annuli 
of ordinary superconductors of thickness d s , separated 
by a layer of oxide of thickness d ox , whose relative di- 
electric constant is e r . Its order-parameter is the relative 
phase angle 4> = Q\ — 62 of the complex order parame- 
ters ^1 = p\exp(i0i) and ^2 = P2<3xp(i02) of the two 
superconductors (labelled 1 and 2). After the transition 
has been implemented, in the adiabatic regime at tem- 
perature T, 4> satisfies the dissipative, one-dimensional 
sine- Gordon (SG) equation 



d 2 (j) 

Th? ' c 2 (T) dt 2 



1 d 2 (f> b d<f> 



c 2 (T) dt X 2 {T) 



sin ( 



(4) 



with periodic boundary conditions ||; x measures the 
distance along the annulus, its width w <C Xj(T) being 
ignored and b is a characteristic frequency that accounts 
for the viscous drag. The velocity c = c(T), which de- 
pends on the nature of the junction, is the Swihart H ve- 
locity, the speed of light in a superconducting-insulating- 
superconducting transmission line. In the Josephson con- 
text, it determines the maximum speed at which the or- 
der parameter <f> can change. 

The topological defects of the JTJ, the solitons of the 
sine-Gordon theory, are termed fluxons. Their static 
equilibrium thickness is the Josephson coherence length 
Aj(T), which plays the role of £, e q{T) earlier. 

Let us attempt to repeat the Kibble- Zurek analysis 
directly on quenching a JTJ with quench time tq. For 
simplicity we begin with a symmetric JTJ, in which the 
electrodes are made of identical materials with common 



critical temperatures T c . At time t after the transition 
Xj(t) = Aj(T(t)) is given by 



Xj(t) = ^h/2ep d e (t)J c (t). 



(5) 



in which J c is the critical Josephson current density. In 
(||) d e (t) is the magnetic thickness. Specifically, if Al(£) 
is the London penetration depth of the two (identical) 
superconducting sheets, then 

d e (t) = d ox + 2Az,(t)tanh 



2A L (t) : 



where 

X L (t) = X L (0)/^/l-(T(t)/T c )* ~ \ L (0)/2^tfc. Nc- 
glecting the barrier thickness d ox <C d s , Xl gives d e = d s 
close to T c , i.e., the magnetic thickness equals the film 
thickness and can be set constant in (|^) . 

All the i-dependence of Xj resides in J c which, for the 
symmetric JTJ has the form 10 



J c (t)^M tanh . A W 



2 ep N 



2k B T{t) ' 



(6) 



In (g) A(t) is the superconducting gap energy and varies 
steeply near T c as 



A(t) ~ 1.8 A(0) 



1/2 




and Pat is JTJ normal resistance per unit area. Intro- 
ducing the dimensionless quantity a — 1.6A(0)/fcsTc 
whose typical valued] is between 3 and 5, enables us to 
write J c (t) as 



J c (t) ~ a,J c (0) (l - = a.J c (0) — 

Thus, in the vicinity of the transition, 

1/2 



corresponding to v = 1/2 in (|l|), where 

6- / * 

V 2e/i O d s aJ c (0) 



(7) 



(8) 



(9) 



On the other hand, for a finite electrode thickness tunnel 
junction, the Swihart velocity takes the form pOj 



c(i) = c yjd ox /t r di{t), 



where 



2 These results, without any additional factors, were origi- 
nally obtained by Zurek on considering the time — t at which 

the field freezes in. 3 A(0) and J c (0) denote the respective values at T = 0. 
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di{t) — d ox + 2Ai(i) coth 



x- 



2X L (t) 



d s \t 



near the transition. Thus c(t) shows critical slowing down 
at the transition, as 



c (*) = c o ( 7- 



1/2 



where Co = co-\/<is^o£c/ e rA|(0). These indices (y = 
1/2,7= 1) are typical of condensed matter systems. The 
causal constraint gives t — a /tqtq, with To = £o/co- In- 
serting reasonable values pj of £0 — 10 fim and Co = 
10 7 m/s, gives To = lps, and assuming tq = 1 s, we find 
i ~ lfj,s. The causal Josephson penetration length is then 

Xj = Xj(t) = £ohf) 2 = £o(^V =10mm. (10) 

This A j, which should characterise fluxon separation 
at a quench for a symmetric JTJ is far too large. For- 
tunately, the manufacture of JTJs typically yields non- 
symmetric devices with more acceptable properties. Sup- 
pose the two superconductors, 1 and 2, now have different 
critical temperatures T C 2 > T c \. Fluxons only appear at 
temperatures T < T c i, from which we measure our time 
t. At this time 



A,(7-,,l- |.sA,(l)l(l-^ 



1/2 



and Ai(i) ~ 1.8 Ai(0) v/ 't / 'tq . The critical Josephson 
current density J' c (t) for a non-symmetric JTJ, being 
proportional to Ai(i)A 2 (t) p0| , behaves just after the 
transition as 



rp \ 1/2 / + \ 1 /' 2 

1 • -€) a M~ Q 



(ii) 



where J' e (0) = ttA 1 (0)A 2 (0)/[A 1 (0) + A 2 (0)]ep N , and 
a' = [Ai(0) + A 2 (Q)]/fcsT c , 1 , provided A 2 (T C)1 ) < 
27rfc_BT c ,i- This is the case here. 

The crucial difference between ( |Tl| ) and ^ is in the crit- 
ical index. Near t = 0, we now find 



Aj(t)=£o 1 



Ik 
T c2 



-1/4 



1/4 



(12) 



where £0 is as in @, since J' c (0) is indistinguishable from 
J e (0) and a' is comparable to a. For the critical behavior 
( |l2| ) to be valid, rather than (||) we need 1 — T c i/T c2 
0(1/ tq) — O(10 -6 ), which is always the case. For a 
typical value (1 — T c \jT C 2) = 0.02 the critical time t is 
now determined by (7 = 3/4) 



t = To /7 T 3 n /7 {l- — 



-1/7 



0.24/xs, 



with our parameters. In turn, 



1.4mm (13) 



is an order of magnitude smaller than Aj of (|io|). 

While new experiments are required, old experiments 
on JTJs by one of us |ll| are compatible with these 
predictions, although their specific parameters are not 
optimal. In these experiments non-symmetric annular 
Nb/Al - AlOx/Nb JTJs (T C)2 /T C) i - 1 « 0.02) with cir- 
cumference C = 0.5 mm were quenched with a quench 
time tq = 0{ls). The intention was, primarily, to pro- 
duce fluxons for further experiments, and the density 
at which they were produced was secondary. From the 
parameters quoted in pl| for sample B, we estimated 



£0 ^ 6.5 /im, c 




s and tq 



0.65ps. Inserting 
these specific values in ( |13| ) gives Aj ~ 1mm (with ex- 
perimental uncertainty of up to 50%). Although C ~ Aj 
we would have expected to see a fluxon a few percent of 
the time, given that the variance An in the number of 
fluxons is A0/27T. Indeed, in practice (invariably single) 
defects formed once every 10-20 times. 

We have no detailed knowledge of how the cooling 
takes place, but do not expect temperature inhomo- 
geneities to be important. The critical slowing down of 
c(t) provides a necessary condition for defects to survive 
inhomogeneity Q and, with empirically comparable £o> 
To, and tq, the situation is no better or worse for JTJs 
than with any other superconducting system undergoing 
a mechanical quench. Other samples of the same circum- 
ference but with different Aj have been tested. Although 
none had C/Xj large, it was observed that the likelihood 
of seeing a fluxon was greater the larger its value, as we 
would have predicted, although this was not quantified. 
This suggests that temperature inhomogeneities are not 
the direct cause of the observed fluxons. 

There are theoretical, as well as experimental, uncer- 
tainties. The SG equation (|J) can only make sense once 
the individual superconductors have adjusted themselves. 
Repeating Zurek's analysis of the Gross-Pitiaevsky equa- 
tion^] for individual superconductors 0] gives a mini- 
mum time at which the sine-Gordon equation is valid 
of is = \/tqTq where To in (||) is now determined 
from Gorkov's equation as To = nh/16kBT c w 0A5ps 
for T c w 10 K. The resulting is ~ OAfxs is commensu- 
rate with the values of i for the typical symmetric and 
non-symmetric cases, falling between them. Whereas this 
suggests that the SG equation is valid for symmetric JTJs 
at time i, it also suggests that we should evaluate \j(t) 
at is, rather than i for the non-symmetric case. However, 
for our typical parameter values the difference between 
A, /(F) and 



4 Justified here by the success of Feynman's coupled model 
equations for ^q and ^2- For example, see Ref.10. 
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\j(t s )^Co 1 



2k 

Tc2 



-1/4 



1/8 



1.1mm (14) 



is so small as to be ignorable, given the crudity of the 
bounds. For the specific sample B of jll| the decrease is 
similar, at Xj(ts) — 0.7mm, and equally ignorable. 

Further, although the prediction together with 
£de/ = has been taken, without additional qualifi- 

cation, as the direct basis for the successful experiments 
HD in 3 He, and experiments in 4 i/e Jl^Jl^l and high- 
ly superconductors Jig ], the causality argument that we 
have presented here is very simplistic. For superfluids 
obeying time-dependent Ginzburg-Landau (TDGL) the- 
ory (and QFT) we know |l6) that, at early times, the 
length £ is, correctly, the correlation length of the fields 
when they have frozen in after the transition. However, 
we also know that the separation of defects is de- 

termined largely by the separation of the zeroes of the 
fields which define their cores. The separation of zeroes 
is a function of the short-range behavior of the correla- 
tion functions |l^,|l8| , rather than the long-range behav- 
ior that determines £. A priori, £ does not characterize 
defect separation. 

Nonetheless, several numerical and analytic cal- 

culations p6| , ^l| , p2[ , based on TDGL theory, have con- 
firmed that the critical index of (||) is, indeed, the cor- 
rect behavior for defect separation. The reason why this 
is so is essentially a matter of dimensional analysis. The 
density of zeroes is, approximately, a ratio of moments 
of the power in the field fluctuations, at early times at 
least. This leads to strong cancellations of the effects of 
the microscopic interactions of the system in question. 

At the same time, the critical time t characteristically 
underestimates the time at which the order parameter 
achieves its equilibrium magnitude, which is a more sen- 
sible time to begin to count defects. However, if these 
other systems are a guide |l^ , ^l| , p2f the true time t*\s 
t* = 0(F), since the unstable long wavelength modes 
that set up large scale ordering have amplitudes that 
grow exponentially. As a result any new scales only oc- 
cur logarithmically in t* jt. In fact, a limited calcula- 
tion, with 4 He in mind, suggests Q that £de/(**) ~ £ 
when counting topological density on an annulus. Thus, 
although the causality bounds are not saturated, their 
consequences (0) and (|j|) survive qualitatively and jus- 
tify experimental confirmation. 

That causality is now seen as a constraint, but not the 
microscopic mechanism, helps explain why the most re- 
cent 4 He experiment |l|] failed to see any vortices. A 
major reason (nothing to do with causality) is that vor- 
tices are most likely to decay much more rapidly p^ ] 
than they were originally thought to do. However, be- 
cause the A He quenches take place entirely within the 
Ginzburg regime, thermal fluctuations make individual 
defects scale dependent [ 21p^ ], and simple dimensional 
analysis most likely breaks down, in a way that causal- 
ity would not have suggested. This is not the case with 
superconductors, for which the Ginzburg regime is very 



small. 

For that and other reasons, we are sufficiently opti- 
mistic to be currently examining the feasability of fab- 
ricating JTJs with larger values of C/Xj with which to 
perform new experiments. 
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